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Abstract. The Brimstone Canary is a geographical highly variable species, distributed throughout southern and eastern 
Africa. Here 1 present data on body measures and plumage characteristics of 476 skins from all parts of the species’ ran¬ 
ge (Fig. 1). In both univariate and multivariate analyses the wing and beak measures showed the most remarked diffe¬ 
rences between populations (Tab. 2, 5). Variation parallels climatic trends with individuals from hot and humid regions 
have shorter wings and smaller beaks (Tab. 7). In contrast to earlier studies plumage colour should not be used to distin¬ 
guish populations, because of high individual variation (Tab. 8). 

Only skins from South Africa (including Natal and Zululand) could be separated clearly (stronger beaks) from all other 
populations (shorter wings, less strong beaks). Clinal variation and high individual variation did not allow more diffe¬ 
rentiation in the northern part of the range. Therefore, the data are in agreement with the separation of the subspecies 
sulphurata (in South Africa) and sharpei (all other populations), while wilsoni and all other subspecies are not suppor¬ 
ted. 

Keywords. Multivariate analysis, subspecies,-wing dimensions, beak size, individual variation. 


1. INTRODUCTION 


The traditional way to describe spatial variation of species 
is the recognition of subspecies. Subspecies are more or 
less clearly separated from other populations of the same 
species by differences in at least one trait and are distrib¬ 
uted in a more or less clearly defined subarea of the 
species' range. This method gives, however, an incomplete 
picture of variation, because many traits vary in dines and 
the spatial characteristics of dines may differ from trait 
to trait (see Gould & Johnston 1972 and Zink & Rem- 
sen 1986 for reviews). 

Hence, to distinguish a representative of a subspecies from 
an individual variant or a local form in a dine profound 
knowledge of individual and geographical variation 
within the species is necessary. Such knowledge was not 
available when most subspecies have been described and 
is still not available for many species. 

The Brimstone Canary Crithagra sulphurata is such a ge¬ 
ographical highly variable species. Within the Canaries 
five afrotropical species form the monophyletic genus 
Crithagra as suggested by morphological, ethological and 
molecular data (van den Elzen 2000). C. sulphurata in¬ 
habits grasslands and savannas with scattered trees 
throughout southern and eastern Africa (Fig.l). Within this 


range and especially in its eastern part a high degree of 
morphological variation has led to the description of sev¬ 
en subspecies: 

-sulphurata (Linnaeus, 1766) (Linnaeus 1766, Syst. Nat., 
ed. 12,1: 305), type from Cape of Good Hope, 

- sharpii (Neumann, 1900) (Neumann 1900, J. Omithol. 
48: 287), t.t. Marangu, Kilimanjaro, 

shelleyi (Neumann, 1903) (Neumann 1903, Omithol. 
Mber. 11: 184), t.t. Kafuro, NW Tanzania, 

-frommi (Kothe, 1911) (Kothe 1911, Omithol. Mber. 19: 
71), t.t. Namanjera, Ufipa, SW Tanzania, 

- loveridgei (van Someren, 1921) (van Someren 1921, 
Bull. Br. Omithol. Club 41: 114), Lumbo, northern 
Mozambique, 

- wilsoni Roberts, 1936 (Roberts 1936, Ann. Transv. 
Mus. 18: 216), t.i. Kloof, Natal, and 

- languens Clancey, 1962 (Clancey 1962, Durban Mus. 
Nov. 6: 193), type from Manhica, Sul do Sava, south¬ 
ern Mozambique. 

Different authors have adopted controversial positions 
which of these types represent real subspecies, e.g.. Fry 
& Keith (2004) in their recent review recognized only 
three of these (sulphurata, sharpei, wilsoni). Rand 
(1968) and Clancey (1972) reviewed available measure¬ 
ments and distribution data and concluded, in accordance, 
that most of the “subspecies” mentioned above are not 
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sharply separated from each other in the traits considered. 
Due to larger body size and stronger beaks South African 
nominate sulphurata can be distinguished from all other 
forms with relative ease. Therefore, Rand accepted only 
two subspecies: sulphurata and sharpei (all populations 
except South Africa). In an alternative view he divided 
sharpei on the basis of subtle differences in plumage 
colour and size into four subspecies: mlsoui (Natal, south¬ 
ern Mocambique), sharpei (Kenya), froiiuui (Tanzania, 
Zambia and Angola) and shelleyi (all other populations). 
In contrast, Clancey (1972) separated four subspecies in 
South Africa alone (sulphurata , wilsoui, kuigueiis and 
shelleyi). Both studies are, however, based on only a few 
traits and from some regions only a small number of in¬ 
dividuals have been included. 


black, of colour, green (negatively correlated to propor¬ 
tion of yellow due to construction of the colour system), 
and red (C. flaviveutris only). Not all traits could be tak¬ 
en for all individuals because of skin condition. Therefore, 
sample sizes differ from analysis to analysis. 

Skins with geographical origin close to each other have 
been pooled for analysis, resulting in 15 ‘populations’ (Fig. 
1). Individuals which locality information were lacking, 
unclear or doubtful (e.g. Collection Meinertzhagen; cf 
Knox 1993) were included in the analyses, but considered 
as belonging to none of the populations. To get larger sam¬ 
ple sizes some of the populations were pooled to larger 
‘population groups’ (A G) in some analyses (Fig. 1). 


In this study 1 present data on a larger variety of meas¬ 
ures and plumage characteristics of C. sulphurata skins 
from all pails of the species’ range. The aim of this study 
is to describe the extent and pattern of geographic varia¬ 
tion in this species. Which of the subspecies described ear¬ 
lier represent real taxonomic units and which are interme¬ 
diate stages arbitrarily picked out of a continuous cline of 
variation? In a second step 1 want to relate morphologi¬ 
cal variation to possible causal factors such as ecological 
conditions (climate, vegetation, topography) and palaeo- 
graphic events. 


2. MATERIALAN1) METHODS 

I measured a total of485 museum skins of Crithagra sul- 
phurata , of which all 476 adults (253 males, 171 females 
and 52 unsexed) were considered in the following analy¬ 
ses. Additionally, I studied 24 skins of the sister species 
(van den Elzen 2000), Crithagra flaviveutris, for compar¬ 
ison (15 males, 9 females). The skins are kept by Zoolo- 
gisches Forschungsmuseum Koenig (ZFMK), Bonn, Zo- 
ologischcs lnstitut und Museum (ZMH). Hamburg, Zool- 
ogisk Museum (ZMUC), Copenhagen, Zoologische 
Staatssammlung (ZSM), Munich, Natural History Muse¬ 
um (BMNH), Bird Group, Tring, and Royal Museum of 
Central Africa (RMCA), Tervuren, respectively. 

Tarsus and beak measures (length, width and height in 
both cases) were taken from all specimens. Length of 
wing, tail and first to ninth primaries as well as Kipp’s dis¬ 
tance (wing tip to first secondary) and graduation of tail 
were measured in all except molting specimens. From the 
wing measures I calculated Wing Pointedness Index I and 

II after Mlikovsky (1978) and Hedenstrom (1989), re¬ 
spectively. Plumage colours (throat, breast, belly and back) 
were compared with standardized colour maps (Scandi¬ 
navian Colour Institute 2001). For statistical analysis I de¬ 
scribed standard colours by four variables: proportion of 


f g 



Fig. 1 . Map of collection localities of skins (crosses - single 
skins; numbers - series with number of collected individuals). 
Ellipses define “populations”. “Population groups” are marked 
by letters A to G. The present breeding distribution of C. sulp¬ 
hurata (simplified after Fry & Keith 2004) is shaded. 


Predominant vegetation structure was revealed for each 
skin locality using the vegetation map of EC-JRC (2003, 
map 2) summarizing all categories under ‘dominant tree 
layer’ as ‘forest’ and all categories under ‘dominant agri¬ 
culture’ and ‘dominant shrub or grass layer' as ‘open habi¬ 
tat’. Climatic data of the closest weather station were tak¬ 
en from Lebedev (1970). Because of the scarcity of weath¬ 
er stations, up to 63 skins gathered on one weather sta¬ 
tion and, in these cases, I used the mean for each trait. The 
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following data were considered: yearly mean of mean, 
maximum, minimum and range of air temperature per day, 
yearly sum of precipitation, daily maximum of precipita¬ 
tion. Altitudinal effects could not be taken into account 
because for most skins altitudinal information were not 
available with sufficient accuracy. 

Non-parametric Factor Analysis were performed using the 
MS Excel Add-In XLSTAT 7.0 (Addinsoft 2003). For all 
other calculations I used the software package SPSS for 
Windows 11.0 (SPSS 2001). 


3. RESULTS 

3.1. Data preparation 

Some variables were transformed to fit a normal distri¬ 
bution: beak height (reciprocal square root transforma¬ 
tion), first primary length (reciprocal), graduation of tail 
(natural logarithmic) and the tarsus measures (cube). All 
other metric variables did not differ significantly from nor¬ 
mal (Kolmogorov-Smimov One-sample-test, p>0.05). 

The plumage colour variables have been subjected to a 
non-parametric Factor Analysis based on Spearman's rank 
correlation. After Varimax rotation the resulting three fac¬ 
tors with eigenvalues larger than I explain 88 % of the 
variance (Tab. 1). The first factor (hereafter called ‘front 
colour’ or ‘Cl’) is highly associated with belly and throat 

Table 1 . Correlations between plumage colour variables and 
factor scores from a Factor Analysis based on a Spearman's rank 
correlation matrix. Correlations > 0.6 printed bold. 


Factor 

Cl 

C2 

C3 

belly: % black 

-0.367 

0.226 

0.027 

belly: % colour 

0.751 

-0.053 

-0.134 

belly: % green 

-0.740 

0.045 

0.175 

breast: % black 

-0.056 

0.248 

0.943 

breast: % colour 

0.336 

-0.223 

-0.747 

breast: % green 

-0.225 

0.010 

0.620 

back: % black 

-0.147 

0.717 

0.332 

back: % colour 

0.243 

-0.951 

-0.188 

back: % green 

-0.141 

-0.307 

0.026 

throat: % black 

-0.567 

0.160 

0.119 

throat: % colour 

0.723 

0.002 

-0.145 

throat: % green 

-0.776 

0.015 

0.248 

eigenvalue 

4.4 

1.6 

1.1 

% variance explained 

38.7 

22.3 

27.8 


colour; higher values represent individuals with more yel¬ 
lowish and less greenish front plumage. C2 (‘back colour’) 
describes back colour, higher values correspond with more 
black and less colour components in the back plumage. 
The colour of the breast is highly correlated with C3 
(‘breast colour’). Individuals with high values of C3 have 
darker and more greenish breasts while low C3 values in¬ 
dicate yellowish breasts. For further analysis I used the 
factor scores of C1-C3 for each skin. 


3.2. Morphological differences between C. siilpliurata 
and C. flaviventris 

At first I looked at differences in morphological charac¬ 
teristics between C. siilpliurata and its sister species, C. 
flaviventris . A Discriminant Function Analysis (DFA) re¬ 
duces the 21 morphological variables to a single canoni¬ 
cal variable with an eigenvalue of 1.13 and a canonical 
correlation of 73 %. Two beak measures and two wing 
measures are included in the analysis (standardized canon¬ 
ical coefficients in parenthesis): beak length (1.533), beak 
height (0.696), 1st primary length (0.463) and distance 
from wing tip to secondaries (0.343). Values for C. fla¬ 
viventris are lower, thus, skins of this species are small¬ 
er in beak and wing dimensions. 95 % of all skins have 
been classified correctly by the DFA. An additional trait 
for differentiation is a red component in the front 
plumage of C. flaviventris: 13 of 21 skins had a slight tinge 
of red, but none of the C. siilpliurata skins. 

3.3. Differentiation of C. siilpliurata populations 

Most traits differ significantly between populations (Tab. 
2). Wing and beak measures show high levels of differ¬ 
entiation (high significance values) in both males and fe¬ 
males. For a detailed analysis I tested for each trait if there 
is a difference between population groups. Again, the most 
pronounced differences are in the wing and beak meas¬ 
ures (Tab. 3). Of 21 population group pairings 16 (males) 
and 11 (females) differ significantly in at least one wing 
measure, respectively, and 12 (males) and 6 (females) in 
at least one beak measure, respectively. The populations 
of Kenya (group G) and South Africa (A) are separated 
by differences in most traits from all other population 
groups, but not from each other. In the central part of the 
species’ range populations vary hardly. There is a tenden¬ 
cy for females to differ in fewer traits than males across 
the same population group pairings. 

Many of the studied traits correlate with each other (e.g. 
wing length and primary length). Intercorrelations may oc¬ 
cur between all variables. Therefore, I performed a Prin¬ 
cipal Component Analysis (PCA) on all 24 original meas- 
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Table 2. Univariate differences in C. sulphurata traits between populations (one-way ANOVA). ^significant differences on the 
5%-level (corrected for 24 tests: P<0.05/24=0.0021). 


males females 


wing measures 

df 

F 

Sig. 

df 

F 

Sig. 

wing length 

13/215 

19.002 

0.000* 

13/146 

10.548 

0.000* 

1st primary length 

13/192 

8.443 

0.000* 

11/127 

4.040 

0.000* 

2nd primary length 

13/192 

7.286 

0.000* 

11/128 

4.660 

0.000* 

3rd primary length 

13/197 

7.083 

0.000* 

11/129 

6.708 

0.000* 

4th primary length 

13/189 

7.733 

0.000* 

11/124 

4.818 

0.000* 

5th primary length 

13/195 

11.383 

0.000* 

11/129 

6.202 

0.000* 

6th primary length 

13/194 

15.328 

0.000* 

11/131 

11.721 

0.000* 

7th primary length 

13/195 

19.416 

0.000* 

11/132 

11.504 

0.000* 

8th primary length 

13/188 

17.909 

0.000* 

11/131 

10.368 

0.000* 

9th primary length 

13/191 

12.721 

0.000* 

11/130 

7.282 

0.000* 

Kipp 

13/213 

2.797 

0.001* 

13/142 

1.433 

0.151 

wing index I 

9/61 

0.534 

0.844 

9/39 

1.242 

0.299 

wing index 11 

tail measures 

9/61 

0.494 

0.873 

9/39 

1.082 

0.397 

length 

13/216 

7.475 

0.000* 

13/145 

5.812 

0.000* 

graduation 

beak measures 

13/214 

1.230 

0.259 

13/143 

1.025 

0.431 

length 

13/217 

9.009 

0.000* 

13/147 

7.826 

0.000* 

height 

12/104 

11.610 

0.000* 

13/73 

9.792 

0.000* 

width 

tarsus measures 

13/218 

20.355 

0.000* 

12/146 

12.610 

0.000* 

length 

13/216 

4.631 

0.000* 

13/147 

5.756 

0.000* 

width 

13/215 

2.576 

0.002 

13/147 

2.681 

0.002 

height 

Plumage colour 

13/215 

2.650 

0.002 

13/147 

3.632 

0.000* 

front colour 

13/218 

3.169 

0.000* 

13/147 

1.487 

0.129 

back colour 

13/218 

4.592 

0.000* 

13/147 

2.047 

0.021 

breast colour 

13/218 

4.067 

0.000* 

13/147 

3.862 

0.000* 


ures (21 morphometric and 3 colour traits). It reduces the 
data set to five uncorrelated factors that together explain 
65 % of the variance (Tab. 4). The correlations between 
the first factor (PCI) and length of wing, tail and all pri¬ 
maries arc all relatively large and positive (1st primary ap¬ 
pears inverse due to transformation). It can be taken as a 
multivariate measure of size; individuals with high PCI 
values are larger in size. The second factor summarizes 
the three beak measures: higher PC2 values indicate in¬ 


dividuals with stronger beaks (beak height transformed as 
1st primary). PC3 and PC4 represent wing shape and 
breast colour, respectively. The fifth factor correlates neg¬ 
atively with front plumage and positively with back 
plumage. Individuals with high PC5 scores are, therefore, 
darker on their back and more greenish on their front, 
while individuals with low PC5 scores are more greenish 
on their back and more yellowish on their front. 
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Table 3. Differentiation of traits (raw measures) between the population groups A to G. Shown arc traits with significant diffe¬ 
rences on the 5%-level (Mann-Whitney-U-test, significance level corrected for 24 tests as in Tab. 2). The diagonal separates ma¬ 
le (right) and female (left) analyses. 

Abbreviations: 1 - length, w - width, h - height; 1-9 first to ninth primary length. B, C, G: proportion of black, colour and green 
of breast, throat, belly and back plumage, respectively. 



A 

B 

C 

D 

E 

F 

G 

A 


wing: 1, Kipp, 
1-9 

beak: 1, h, w 
tail 1 

breast: BC 

wing: 1, Kipp, 
1-9 

beak: 1, h, w 
tail 1 

wing: 1,2, 8 

beak: 1, w 

wing: 1, 1,3, 6 

beak: 1, h, w 
tail 1 

wing: 1, Kipp, 
1-9 

tarsus: w, h 
beak: 1, w 
tail 1 

belly: CG 
breast: BC 
throat: G 

beak: h, w 

B 

wing: 1 
tarsus: 1 
beak: 1 


wing: 1, 1-9 

wing: 2, 4, 5 

wing: 1, 1-9 
tarsus: 1 

tail 1 

back: G 

wing: 1, 1-9 

beak: w 

throat: G 

wing: 1, 1 -9 
tarsus: 1 
beak: 1, w 
tail 1 

back: G 

C 

Schn: L 






wing: 1, 1-9 

Schn: L , H , 

Br 

D 

wing: 1 

beak:1 
tail 1 




wing: 8 


wing: 1, 2, 7, 8 

beak:1 

E 

beak: 1, h, w 

wing: 1, Kipp, 

2, 4-9 


w ing: 1, 4-9 


wing: 1, 2, 5-9 
tarsus: 1 

tail 1 
belly: G 
breast: C 

back: BC 

wing: 1, 1,3-9 

beak: 1, h, w 

F 

wing: 1 
tarsus: 1 
beak: 1, h, w 
tail 1 

breast: BC 

wing: 7 



wing: 1, 1-9 

tail 1 

breast: BC 
back. BC 


wing: 1, 1-9 
tarsus: 1, w 
beak: 1, h, w 
tail 1 

G 


wing: 1, Kipp, 

2, 7-9 
tarsus: 1 
beak:1 


wing: 1 

wing: 7 

wing. 1, 2, 3, 
6-9 
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Table 4. Correlations between the 24 original measures and principal component scores from a PCA (using Varimax rotation with 
Kaiser normalization) of the correlation matrix of 476 C. sulphurata skins. Correlations >0.600 are printed bold. 

Variable 

PCI 

PC2 

PC3 

PC4 

PC5 

wing length 

0.864 

0.331 

0.137 

0.105 

0.013 

1st primary length 

-0.726 

-0.171 

0.199 

-0.108 

-0.156 

2nd primary length 

0.780 

0.124 

-0.274 

0.087 

0.173 

3rd primary length 

0.800 

0.120 

-0.292 

0.029 

0.120 

4th primary length 

0.822 

0.077 

-0.298 

-0.013 

0.059 

5th primary length 

0.858 

0.119 

-0.236 

0.039 

0.011 

6th primary length 

0.891 

0.208 

-0.012 

0.103 

0.022 

7th primary length 

0.919 

0.215 

0.091 

0.100 

0.040 

8th primary length 

0.905 

0.224 

0.141 

0.092 

0.050 

9th primary length 

0.844 

0.203 

0.137 

0.100 

0.102 

Kipp distance 

0.366 

0.198 

0.234 

0.036 

-0.208 

wing index I 

- 0.111 

-0.029 

0.946 

0.037 

0.041 

wing index 11 

-0.110 

-0.012 

0.943 

0.038 

0.007 

tail length 

0.636 

0.273 

0.081 

0.125 

0.105 

tail graduation 

0.012 

0.182 

-0.035 

0.461 

-0.082 

tarsus length 

0.286 

0.048 

0.052 

0.625 

0.189 

tarsus width 

0.094 

0.369 

-0.015 

0.106 

0.296 

tarsus height 

0.181 

0.294 

0.004 

0.209 

0.319 

beak length 

0.400 

0.697 

-0.010 

0.148 

-0.009 

beak height 

-0.218 

-0.823 

0.033 

-0.020 

-0.023 

beak width 

0.382 

0.729 

-0.004 

0.078 

0.096 

front colour 

0.032 

-0.153 

0.027 

0.133 

-0.661 

back colour 

0.170 

-0.055 

0.050 

0.009 

0.683 

breast colour 

0.036 

-0.011 

0.071 

0.810 

-0.081 

eigenvalues 

8.172 

2.463 

2.270 

1.452 

1.296 

% variance explained 

34.0 

10.3 

9.5 

6.1 

5.4 


Table 5. Standardized Canonical Discriminant Function 
skins from the population groups A to G. 

Coefficients from a 

DFA of the five PCA scores of 443 C. 

sulphurata 

PCA-factor 

DF1 

DF2 

DF3 

DF4 

PCI 

0.824 

-0.633 

-0.272 

-0.007 

PC2 

0.816 

0.612 

-0.32 

-0.184 

PC4 

0.52 

-0.026 

0.779 

-0.455 

PC5 

0.516 

0.178 

0.365 

0.809 

eigenvalues 

1.517 

0.244 

0.136 

0.015 

% variance explained 

79.3 

12.7 

7.1 

0.8 

canonical correlation 

0.776 

0.443 

0.346 

0.122 
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Table 6. Proportion (%) of skins classified correctly as belon¬ 
ging to population group A-G by the DFA. N - sample size. 



N 

A 

B 

c 

D 

E 

F 

G 

A 

29 

79 

0 

0 

0 

3 

3 

14 

B 

27 

0 

0 

0 

0 

0 

100 

0 

C 

36 

0 

0 

0 

0 

8 

89 

3 

D 

14 

0 

0 

0 

0 

21 

79 

0 

E 

72 

3 

0 

0 

0 

40 

53 

4 

F 

231 

1 

0 

1 

0 

3 

94 

1 

G 

34 

9 

0 

0 

0 

18 

21 

53 


To determine which of these factors have the highest se¬ 
lectivity for the population groups, I performed a stepwise 
Discriminant Function Analysis (DFA) on the PCA scores 
(Tab. 5). PC3 (wing shape) does not contribute to differ¬ 
entiation between groups and is removed by the analysis. 
The other four PCA scores are transformed to four DFA 
factors. 65 % of skins are classified correctly by this analy¬ 
sis. There are, however, marked differences between the 
population groups in the proportion of correct classifica¬ 
tion (Tab. 6). Population group A (South Africa) is clear¬ 
ly divided off from all other groups. Population group G 
(Kenya) is separated less, but more than 50 % of skins are 
classified correctly. Most overlapping of G-skins is with 
group A. Nearly all other skins from groups B, C, I), E 
and F are classified as belonging to population group E 
(Congo, SW Tanzania) or F (Kivu/Uganda/Rwanda). Sep¬ 
arating males and females in a single DFA each gives sim¬ 
ilar results and is, therefore, not presented here. 



C 




F 



-3 0 3 


ro 

o 

_Q 

CM 

O 

CL 


PCI (size) 


Fig. 2. Crithagra sulphurata skins in the space defined by principal components 1 and 2 (see Table 4). Solid symbols - all skins; 
filled symbols - skins from population group A to G, respectively. 
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The groups B to F of eastern Africa which appear homog¬ 
enous in the above population group analysis may be dif¬ 
ferentiated on a lower level. Therefore, I calculated a DFA 
at population level on all skins except from population 
group A and G. Again, most skins are classified as belong¬ 
ing to populations of group E or F. 

3.4. Pattern of geographic variation 

The large range of classification accuracy by DFA (Tab. 
6) indicates that some population groups are more sepa¬ 
rated from their neighbours than others. Scores of all 
groups overlap to some extend with the characteristics of 
other populations (Fig. 2). The pattern of variation is much 
more complex than described by clear-cut population 
groups. As an example, in Fig. 3 / Fig. 4 wing length and 


20°W 30°W 40A/V 





82 0-84 9 mm 
0 79 0-81.9 mm 

% 76 0-78 9 mm 
• 73 0-75.9 mm 
* <72.9 mm 

Fig. 3. Variation of wing length of C. sulphurata. Shown are 
means for 2°30'-fields with at least one studied skin. The solid 
line schematizes the coastline of southern and eastern Africa. 


beak width are shown against the geographical origin of 
the skins. 

The population group with lowest overlap to other groups 
is South Africa (including Natal; group A). South African 
C. sulphurata are characterized by strong beaks and large 
to intermediate body sizes (Fig. 2). Within population 
group A both wing and beak measures are relatively ho¬ 
mogeneous (Fig. 3, Fig. 4). Skins from Natal have a lit¬ 
tle bit lower scores in both principal components, thus dis¬ 
tinguished by little weaker beaks and little shorter 
wings/tails. The differences to the next population in the 
north (group C) are, however, much more remarkable than 
the differences between the South African and Natal pop¬ 
ulation (Figs 2-4). 


20°W 30°W 40°W 



8.4-9 0 mm 
7.7-8.3 mm 

0 7.0-7.6 mm 
• 6.3-6.9 mm 
* 5.6-6.2 mm 

Fig. 4. Variation of beak length of C. sulphurata. See Figure 
3 for explanation. 



Bonner zoologischc Beitragc 55 (2006) 


69 


In relation to all other groups Kenyan C. sulphurata (group 
G) look similar to skins from the other end of the species 1 
range, group A. Beaks of skins from Kenya are interme¬ 
diate between South African birds and all other groups, 
while body size is larger in Kenya than in all other groups 
(Fig. 2). Overlap to neighbouring populations is, howev¬ 
er, considerable. As shown by Fig. 3 and Fig. 4 the pop¬ 
ulation of Kenya is less homogenous in both beak and 
wing measures. Skins from northern Kenya have as long 
wings as skins from southern Kenya, but in respect to beak 
size they resemble birds from the eastern and central part 
of the species 1 range. 

Outside Kenya and South Africa no remarkable morpho¬ 
metric differentiation between population groups exist. 
There is a tendency, however, for skins from Mocambique 
(B) and Angola (D) to be smallest in both principal com¬ 
ponents (Fig. 2). Wing length increases from the coasts 
to central Africa and could be described as a quadratic 
function of geographic longitude (WING=-1.23E- 
10*LONG 2 + 6.87E-05*LONG + 66.25; R 2 =0.06; ANO- 
VA df=2/374 F=11.1 P=0.000). Beak width variation is cli- 
nal, too, and fits a linear function with highest values at 
the eastern coast (BEAK=-1.51 E-06*LONG+6.82; 
R 2 =0.03; ANOVA dl—1/376 F=11.2 P=0.001). Apart from 
this clear trend, for both measures areas with low scores 
alternate with high score areas, however (Fig. 3, Fig. 4). 

The clinal longitudinal trend is visible in the PC1/PC2 plot 
as a shift along the axis B (Mozambique) to C 
(Malawi/Zimbabwe) to E (Congo/Tanzania; Fig. 2). Due 
to high intra-group variability both population group F 
(Kivu/Rwanda/Uganda) and D (Angola) overlap with 
nearly all other central population groups. Skins from 
these two groups are, therefore, indistinguishable from the 
other population groups except A and G. Hence, the DFA 
could not dissolve the population groups B to F (Tab. 6) 
and classified all skins from the central part of the species’ 
range as belonging to group F. 

3.5. Impact of climate and vegetation 

To reveal associations between variation of morphologi¬ 
cal traits of C. sulphurata and climatic conditions 1 per¬ 
formed a PCA with the six climatic variables mentioned 
above and the five principal components of morphologi¬ 
cal measures (cf Tab. 4). The resulting two factors explain 
45 % of total variance (Tab. 7). The first principal com¬ 
ponent (cPC 1) shows that body size (PCI) decreases along 
a gradient of increasing temperature. Large C. sulphura¬ 
ta can be found under relatively moderate temperature 
conditions. Beak dimensions decrease as sum of precip¬ 
itation increase. Therefore, strong beaks are associated 
with relatively dry climate. 


Table 7. Correlations between the 6 climatic variables, the 5 
principal components (cf. Table 4; mean for each weather sta¬ 
tion) and principal component scores from a PCA (using Vari- 
max rotation with Kaiser normalization) of the correlation ma¬ 
trix of 68 weather stations. Correlations >0.6 are printed bold. 


variable 

cPCl 

cPC2 

mean of temperature 

0.898 

-0.240 

maximum of temperature 

0.850 

-0.146 

minimum of temperature 

0.916 

-0.089 

range of temperature 

-0.262 

-0.057 

sum of precipitation 

-0.041 

-0.801 

maximum of precipitation 

-0.468 

-0.064 

PCI 

-0.641 

0.443 

PC2 

0.149 

0.785 

PC3 

-0.156 

0.016 

PC4 

-0.089 

0.327 

PC5 

-0.238 

0.329 

eigenvalues 

3.4 

1.5 

% variance explained 

31.0 

13.9 

Table 8. Mean variance of principal components (see Tab. 4) 
at eight localities (“local”; between 8 and 15 skins, total N=59) 
and in a random sample of the respective populations (“popu¬ 
lation”; same sample sizes). “Significant”: number of the eight 
local-population pairings with significant differences in varian¬ 
ce (Levenc’s test, p<0.05). 

Component Local 

Population 

Significant 

PCI (body size) 0.59 

0.93 

2 

PC2 (beak size) 0.48 

0.79 

1 

PC3 (wing shape) 1.03 

0.53 

1 

PC4 (breast plumage) 0.77 

1.19 

1 

PC5 (front/back plumage) 0.80 

1.44 

1 

Skins from forested areas have lower scores of the size- 


related first principal component (thus, are smaller in size; 
mean -0.12) than skins from open landscape (mean 0.07; 
t-test t=2.06 clf=A\6 P=0.040; P all other PC >0.05). This 
difference is more pronounced in Kenyan skins (popula¬ 
tion group G), where it is highly significant (t=3.76 df= 32 
P=0.001). 
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3.6. Individual variation 

From seven localities I have series of eight or more skins 
that had been sampled in the same season (within three 
month). I compared variability of these samples with vari¬ 
ability in a random sample of same sample size taken from 
the same population. In principle variances of local sam¬ 
ples are slightly smaller than population variances (Tab. 
8), but this difference is significant in only one or two of 
eight pairings. Mean variances of the eight samples are 
smaller in principal components 1 and 2 (body and beak 
size) than in the other variables (plumage colour). This is 
true for both local and population level samples. 

4. DISCUSSION 

The morphological variability of C. sulphuratci is consid¬ 
erable. Differences in body dimensions and plumage 
colour occur not only between populations but also indi¬ 
vidually at a single locality (Tab. 8). Therefore, if geo¬ 
graphical variation is studied attention has to be paid to 
individual variation, too. 

In both univariate and multivariate analyses wing and beak 
measures show the most striking differences between pop¬ 
ulations (Tab. 2, 5). Wing and beak measures also exhib¬ 
it lowest individual variation at one locality (Tab. 8). 
Therefore, these traits are suitable to separate populations 
and subspecies. The revisions of Rand (1968) and 
Clancey (1972) were based on wing and beak measures, 
too. Additionally, both authors considered slight differ¬ 
ences in plumage colouration (a little more greenish/ycl- 
lowish etc.) to recognize subspecies. My study of a larg¬ 
er sample of skins shows, however, that plumage colour 
varies enormously at one site (Tab. 8) and can not be used 
to distinguish populations. 

Morphological variation of C. sulphurcita parallels climat¬ 
ic trends: individuals from hot and humid regions have 
shorter wings and smaller beaks (Tab. 7). The pattern of 
variation in wing and beak size (Fig. 3, Fig. 4) could be 
explained by altitudinal effects to a large extend. Similar 
climatic conditions may be responsible for similar trait di¬ 
mensions of the populations of Kenya and South Africa, 
resident at the opposite borders of the species' range. Dif¬ 
ferences in body size in relation to predominate vegeta¬ 
tion cover are probably a side-effect of climate, too, be¬ 
cause also in rain forest areas C. sulphurata settles in the 
scattered open habitats and not in forests (Fry Sc Keith 
2004). Especially in Kenya, where C. sulphurata inhab¬ 
its a broad range of altitudes, the relation between body 
size on one hand and climate and vegetation on the oth¬ 
er hand is strong. In this population separation of sites 
close to each other but differing in altitude is obvious. If 


wing length is taken as a measure of body size, this cli¬ 
matic trend parallels Bergmann’s rule, which suggest an 
increase of body size from cold to warm climates (Zink 
Sc Remsen 1986). 

There is little evidence for geographical isolation of pop¬ 
ulations due to climatic changes in the past 20.000 years. 
Palaeovegetation records certify distinct changes in veg¬ 
etation cover during this period (Adams 2004). Habitats 
occupied by C. sulphurata today (grasslands, scrub, sa¬ 
vannas) have been driven back due to moister conditions 
during the holocene (about 10.000 years ago), resulting 
probably in a much smaller range of C. sulphurata at that 
time. However, open landscapes ranged continuously from 
southern Africa to Kenya, at the most interrupted by a nar¬ 
row forest belt at 20-25° S. 

Geographic variation reflects adaptation to different en¬ 
vironmental conditions within a species’ range. Therefore, 
it could be seen as a model of evolutionary steps of a sin¬ 
gle population over time in a changing environment 
(Gould Sc Johnston 1972). Differences in wing shape and 
size between populations often parallels migratory behav¬ 
iour (Leisler Sc Winkler 1985). Fry Sc Keith (2004) de¬ 
scribe C. sulphurata as resident throughout its range. Short 
distance (altitudinal?) movements may occur in some pop¬ 
ulations, but at the most occasionally (Mackworth-Pread 
& Grant 1963, Harrison et al. 1997). Therefore, mor¬ 
phological differences more probably reflect differences 
in habitat structure and use since body size and shape af¬ 
fect manoeuvrability (Zink Sc Remsen 1986). As Leisler 
& Winkler (1985) have shown in their comparative study, 
morphological characters are strongly correlated with for¬ 
aging techniques. The variability of beak dimensions in 
this study suggests differences in foraging strategies or 
food composition of habitats within the range of C. sul¬ 
phur ata. In all parts of the range C. sulphurata settles in 
open habitats with scattered trees, food mainly consists of 
seeds and small fruits (Mackworth-Pread Sc Grant 
1963, Harrison ct al. 1997, Fry Sc Keith 2004). How¬ 
ever, detailed information on geographic variation of habi¬ 
tat and food preferences is lacking. More field studies of 
the species are necessary to reveal this relationship. 

Single traits can be linked due to developmental con¬ 
straints. The PC A shows that the different measures of 
wing and tail size do not differ independently, but repre¬ 
sent a single principal component (Tab. 4). This corre¬ 
sponds with many other studies, where the first, highly 
correlated factor could be interpreted as a size factor 
(Gould & Johnston 1972). The fact, that beak dimen¬ 
sions, wing shape and plumage colour each form separate 
factors, points to habitat or food related adaptations inde¬ 
pendent of size. 
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Only the skins from South Africa can be separated clear¬ 
ly from all other populations with the data presented here. 
The clinal nature of variation, the large overlap of trait 
measures and the large individual variation do not allow 
the recognition of subspecies in the remaining, northern 
part of the species' range. 

C. sulphurata of the South African nominate subspecies 
sulphurata are characterized mainly by their strong beaks 
(Fig. 2). Tlie range of this subspecies includes South Africa 
from the Cape Province to Natal and Zululand. In contrast, 
Rand (1968) recognized the population of the Cape 
Province as a separate subspecies, too, but classified the 
populations of north-eastern South Africa (Natal/Zululand) 
as belonging to another subspecies which includes all oth¬ 
er populations in the north. Already Mees (1970) point¬ 
ed to the fact, that the populations of the Cape Province 
and Natal are different from each other in plumage colour 
at the most, but not in other measures (see also Rand's 
own measurements). He therefore argued for combining 
both populations into one subspecies. The statistical analy¬ 
sis of skins presented here allows no separation between 
individuals from the Cape Province and Natal/Zululand, 
but separates between these and the populations in the 
north. These are isolated geographically, too, because C. 
sulphurata is lacking along the Limpopo River (Harri¬ 
son et al. 1997). 

Individuals distributed in the northern part (sharpei) are 
- despite marked individual variation - smaller (shorter 
wings) and have less strong beaks than their South African 
congeners. Both traits vary outside South Africa in a com¬ 
plicated pattern without abrupt changes. Variation in both 
wing length and beak size have a clinal and an altitudi¬ 
nal compound. The separation of subspecies is arbitrari¬ 
ly and would not contribute to a better understanding of 
geographical variation in this case. At best the population 
of the Kenyan highlands could be recognized by relative¬ 
ly strong beaks, but the transition to neighbouring popu¬ 
lations is fluid. A possible differentiation in traits not stud¬ 
ied here (e.g., genetical or ethological) have to be reserved 
for future studies. 

In conclusion, the data presented here are in agreement 
with the viewpoint of Fry & Keith (2004) in respect to 
the subspecies sulphurata and sharpei , while the recog¬ 
nition of a subspecies wilsoni is not supported by this da¬ 
ta set. A further differentiation of populations in the north¬ 
ern part of the range as done by other authors (see above) 
could not be confirmed in this study. 
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